arXiv: 1502.02738vl [math.PR] 9 Feb 2015 


On the range of the transient frog model on Z 


Arka Ghosh* Steven Noren^ Alexander Roitershtein* 


February 11, 2015 


Abstract 

In this paper, we observe the frog model, an infinite system of interacting random 
walks, on Z with an asymmetric underlying random walk. Under the assumption 
of transience with a fixed frog distribution, we construct an explicit formula for the 
moments of the lower bound of the model’s long-run range, as well as their asymptotic 
limits as the drift of the underlying random walk vanishes. We also provide conditions 
in which the lower bound can be scaled to converge in probability to the degenerate 
distribution at 1 as the drift vanishes. 
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1 Introduction 

Consider the following interacting random walks model on Z: initially at each site x there 
is a fixed number rj x of sleeping particles (“frogs”), and there is a certain number (r/ 0 ) of 
active frogs at the origin. The active frogs perform in discrete time, simultaneously and 
independently of each other, a biased (say, to the right) nearest-neighbor random walk on 
Z. When an active frog visits a site x, it activates the r] x sleeping frogs at x, in which each 
active frog performs the same underlying random walk starting from its initial location, all 
random walk transitions being independent of each other. The active frogs continue to visit 
other sleeping frogs and activate them. This model for an infinite number of interacting 
random walkers is called the frog model on Z (with drift). 

Following ra. we define the notions of recurrence and transience for the frog model as 
follows: 

Definition 1.1. A frog model is called recurrent if 0 is visited infinitely often by active frogs 
w.p.l, and transient if 0 is visited only finitely often w.p.l. 
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It is shown in [12] that the zero-one dichotomy actually takes place, namely a one¬ 
dimensional frog model is either transient or recurrent. Both necessary and sufficient condi¬ 
tions for recurrence of the frog model on Z based on the configuration of frogs and the drift 
of the random walk are provided in [T2J. Recurrence for variants of the frog model on more 
general graphs have been first explored in [26] (for the symmetric random walk on 7L d ) and 
subsequently in 0. m and [23]. Shape theorems for the model in 7L d have been obtained 
in m Ej. For further background on the frog model and its variants, refer to [23]. For an 
account of the most recent activity in the area see E ei nans]. In particular, [ 8 ] generalizes 
a recurrence criterion of [T2] to a model in Z rf , [T5] and [T 6 ] provide recurrence and transience 
criteria for the frog model on trees, and |5] studies survival of particles in a one-dimensional 
variation of the model, also partially extending some of the results of [ 12 ] , 

The frog model can be interpreted as an information spreading network j2]|23]. The 
underlying idea is that an active frog holds some information and shares it with sleeping 
frogs when they meet, activating the sleeping frogs who then spread the information along 
their random walk path. A closely related to our model particle process on Z, describing the 
evolution of a virus in an infinite population (e.g., computer network), has been considered 
in pL 8 l [T9] and [5]. The model is also a discrete-time relative of the one-dimensional stochastic 
combustion process studied in [7] [23] - 

In this article, we will explore the behavior of the frog model, in particular its range, 
when transience is assumed. We specifically introduce a drift component to the random 
walk and explore how its magnitude effects the range of visited sites in the model. Each 
active frog will move one integer to the right with probability p G (|, 1 ), or one integer to 
the left with probability 1 — p. Thus the underlying random walk is transient to the right. 
We define the drift constant 

p:= —e(0,i). (1) 

V 

The drift term p can be seen as a measure of “transience” of this frog model; small values 
of p indicate more frequent rightward movement by the frogs, whereas values of p close to 1 
more closely resemble recurrence with a slight rightward drift (see, for instance, formula (J 8 ]) 
below for a concrete random walk result). Of particular interest is the collective behavior of 
the frogs as p t" 1 . 

By Theorem 2.1 in [12], this frog model is transient when we assume an identical distribu¬ 
tion of frogs on the nonnegative sites. It is easy to see that w.p.l there must be only a finite 
number of visited sites to the left of the origin. That is, in the language of [3], transience 
implies local extinction for our model. 

In this article, while assuming r) x — 0 for all x < 0, we will first explore the single-frog 
case, i.e., the frog model in which rj x — 1 for all x > 0 (cf. [5, IS[19]). We will provide exact 
and asymptotic results for the moments of the lower bound of the range, which will be used in 
convergence theorems. After that, we will move to more general choices for 77 and show that, 
under certain conditions, the frog model’s lower bound will behave asymptotically similar to 
that of the single-frog case. Finally, we will provide asymptotic bounds for moments of the 
frog model range when 77 is supported on all of Z. 

The rest of the paper is organized as follows. A short Section [2] introduces notations and 
certain technical tools necessary for our proofs. The three subsequent sections constitute 
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the main body of the article. The single-frog case is considered in Section |3j A class of 
more general initial configurations of frogs g is discussed in Section [4j The consideration of 
configurations supported on Z is discussed in Section 0 

2 Preliminaries 

For our calculations of the density and moments of the random variable representing the 
range of the frog model, we make use of common notation in analytic number theory and 
combinatorics: 

Definition 2.1. For all ci,q such that \q\ < 1, the q-Pochhammer symbol is defined by 


c—1 

(a; q) c := ]^[(1 — aq j ) for c G N U { 00 }. 
j =0 

The g-Pochhammer symbol is one of the key functions in the construction of g-analogs 
in number theory, and is often used in the theory of basic hypergeometric functions and 
analytic combinatorics muni Eg. For example, JT3] provides the following pair of power 
series identities: 


( 2 ; 9)oo 


OO 


E 

n =0 


)rign(n-l)/2 

( 9 ; q) n 


X 


( 2 ) 


1 

(259)00 


OO 


E 

n =0 



(3) 


The power series above reveal that, for fixed nonzero q E (—1,1), g(z) := (z;q) 00 and 
h(z) : = are both analytic functions on (0,1). Hence, for any analytic function 

/ : (0, 00 ) —» (0,1), the compositions g o / and ho f are both analytic functions, and 
thus infinitely differentiable. 

For the calculation of generating functions, it is also useful to recall the following q- 
Pochhammer equality: 

(q x+1 -, 9)00 = X EN (4) 

( 9 ; q)x 

Definition 2.2. Let q E (—1,1). The q-gamma function is defined as 


r ff (*) 


( 9 ; 9 )OO 

(9 z ;9)oo 


(i - 9) 1 " 2 


The q-digamma function ip p is defined as 


M z ) 


1 dT q (z) 
r g(z) dz 


OO 

= — ln(l — q) + In q 

n =0 


~n-\-z 


-fTl+Z 


(5) 
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To facilitate our calculation of the range moments, we need to develop notation for Bell 
polynomials. The Bell polynomials are defined (see, for instance, mm) as the triangular 
array of polynomials m > k, given by 

B m ,k(X I; *^2) • • • j •En—k+ 1) 

\- ml pi\ fcl p 2 \ fca / a^m-fc+i \ fcm ~ fc+1 

ki\k 2 \ ■ ■ ■ k m - k+ i\ V 1! / V2!/ \(m — k + 1)\) 

where the sum is taken over all sequences {fc 1; k 2 , ■ ■ ■, k n -k+ 1 } C N U {0} satisfying 

h + k 2 H-1- k m - k+ 1 = k and k ± + 2k 2 H-h (m - k + l)k m _ k+1 = m. 

We define the m th complete Bell polynomial to be 

m 

B m (xi,x 2 , ..., x m ) = E Bm,k (xi,x 2 , • • • , 1) • 

k=1 

It will be important for our asymptotic analysis of the range of the model to note that 
B m (x i, X 2 , • • •, Xm) has only one term of the form x{, namely x™, and that the coefficient for 
this term in B m (x i, x 2 ,..., x m ) is 1. 

While the Bell polynomials has many intriguing details that can be explored in combina¬ 
torial number theory (see, for instance, [6] and references therein) , we will mostly concern 
ourselves with their presence in the celebrated Faa di Bruno’s formula [13 ng for higher 
derivatives of composite functions: 

( 7 \ rn rn 

j) = B m , h (g\t),g''{t),...,g^- k+1 \t)). (7) 

' k= 1 

Throughout this paper, we make use of the common asymptotic notation: f(x) ~ g{x) as 
x —y c if and only if lim x _>. c = 1. 

3 The single frog per site case 

We will first assume that r/ x = 1 for each nonnegative integer x and tj x — 0 elsewhere. Let 
p E (0,1) and consider the corresponding frog model. Let W p represent the random variable 
for the negative of the minimum of the visited sites in this model. For convenience, we will 
construct a family of mutually independent random variables (A" p ) pe ( 0 ,i) that all share the 
same probability space (hi, J 7 , P ) such that X p and W p share the same distribution. 

The distribution function of X p can be easily found by observing that, by the rightward 
tendency of the initial active frog, all frogs w.p.l will eventually be woken. Furthermore, if 
(*S'n)n>o has the distribution of the underlying random walk, then 

P(S n = 0 for some n > 11 S 0 — 1) = p, (8) 

where p is introduced in ([T]). With this observation, we see that for all x > 0, 

OO 

P(x p < x) = P(X p < X +1) = n(i - P x+J ) = (P* +1 ; p) oo- (9) 

3 = 1 
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It is a simple exercise to see that P(X p = 0) = J^T^l — P*) — (P'i P) oo- For all x > 0, the 
value of Xp s probability density function is 

OO 

P(X p = x) = P(X p <x)~ P(X p < X - 1) = J](l - p x+j ) 

3 = 1 
oo 

= (i - (i - p*» na - =p’(p I+ t pu - 

3 = 1 

With the density known, we would now like to study the behavior of X p for values of p close 
to 1, where the frog model more closely resembles the recurrent case. Objects to observe the 
concentration of the distribution of X p include the central statistics of the random variable, 
such as the mode and the expectation. For fixed p, we define a unique representative of the 
mode statistic, named M p , by 

M p := min{x > 0 : P{X p — x) > P(X p = n) for all n > 0}. 


IF 1 - p x+i ) 

3=0 


( 10 ) 


While the mode statistic is not usually observed compared to other central statistics of a class 
of random variables, it is a quick calculation that can often provide insight on asymptotic. 
Also, as we observe later, the concentration of X p around its mode will be influential to its 
limiting behavior. With that in mind, we present the following result: 


Theorem 3.1. Let p G (0,1), and let M p be the mode statistic for the distribution of X p . 
Then, 


ln(l — p) — In p 
In p 


<M p < 


ln(l — p) — ln(2 — p) 
In p 


where |_uj is the largest integer less than or equal to a G M. In particular, M p ~ ln ^ ^ as 

p^l. 


Note that for all p G (0,1), the difference between the two bounds in the theorem’s 
conclusion always belongs to the open interval (0,2). Hence, even for values of p that make 
the bounds extremely large, the theorem narrows down M p to two possibilities. 

Proof of Theorem \3 . 1\ Considering x to be a continuous variable, we note that 


d_ 

dx 


(p x+ \ 


P) 


d 

oo = 3 - exp 

dx 


OO ,00 

(£ Mi - p r+i )) = ( E Mi - p i+j )) (p I+1 ; p) 


(ii) 


assuming that the derivative of the series exists. By Theorem 7.17 in [25], the existence 
of the derivative depends on the uniform convergence of the series of derivatives on closed 
intervals. It is enough to confirm that Y1JL] 11 i-^+j I loo absolutely converges for any closed 
interval [a, b] we choose, where ||/||oo = su P{|/( a; )l : x £ The series is bounded 

above by 2 1 so it converges absolutely. Hence, the series in (HI]) is differentiable, and 
furthermore, by [25], Theorem 7.17], the derivative of the series is the series of the derivatives. 
Thus, 

<7 EE n x +3 

^(p x+ 1 ; p )00 = - Hp)(p x+1 ; p)oo 

3 =1 ^ 
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Using the above derivative, we can find the derivative of the density function for X p : 


d_ 

dx 


x+1. 



sUlo 

3 = 1 


P X+J )P 2 


ln(p)p X (p X+1 ,p)oo(l - 

3 = 1 


P X+J \ 

1 - p- x +J / ‘ 


For the given p, ln(p)p x (p x+1 , p)^ is nonzero for all x > 0. Therefore, by the l st -derivative 
test for critical points, M p is an integer within 1 away from the positive value m such that 


OO 


£ 


pm+j 
1 _ p m+j 


1. 


Note that p x+J < ^ t+i f° r all x > 0 and all j E N. Thus, using the formula for 

the geometric sum, 


P 


.m+1 


p 


£p” +J <£ 


P 


m+j 


P 


i m+j 


- = 1 < V 
1 ~ 


P 


m+j 


P 


m +1 


j =1 j =1 j =1 

The result follows accounting for the fact that M p is integer-valued 


l-p m+1 l-p m+1 l-p‘ 


□ 


We would now like to compare the mode of X p with the moments of the random vari¬ 
able. Calculating the moments directly from definition can be quite a challenge, given the 
convoluted expression of the density in (OTT) . Therefore, we take the more circuitous option 
of first finding the cumulants of X p . 

We first define the cumulant generating function of X p . With the moment generating 
function M Xp (t) E(e tXp ), the cumulant generating function g p (t) of X p is defined as the 
natural logarithm of 

9pit ) := log (M Xp (t)) = log (E(e tXp )). 

We then define the m th cumulant k^ of X p to be the m th derivative of the cumulant 
generating function evaluated at 0: 



Cumulants are known in probability to be an alternative to moments which are difficult to 
calculate explicitly. Cumulants can also be used to determine moments through the use 
of Faa di Bruno’s formula (J7J) (see, for instance, [201, [21] and references therein), which 
is the direction we will take. It turns out that the cumulants of X p , though unable to 
be written down using fundamental functions, can be expressed as straight-forward series 
representations. 

Lemma 3.2. For each p E (0,1), the cumulant generating function of X p is 


9 P {t) 


k =1 


l-p fc 


e tpk 


Furthermore, for each m E N, the m th -cumulant of X p is 

°° Um —1 


- 

k p ~ 


k=\ 


P‘ 


k 


( 12 ) 


(13) 
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Proof. Using (|3l) . we calculate the moment generating function of X p : 


OO tv v 

e p 


(p; p) c 




i-P fc 

L± 1 - eV fc ' 


(14) 


Taking the natural logarithm of gives us 
derivative of g p (t): 


To find we first find the m 


th 


oo ^7 

»?■><*> = E(* 

fe=i 
00 00 


In 


1 ~p k 
1 - e t p k 


OO 

= £' 

fc=i 

OO OO 


/d 


EU) (mi-^)+E 


(eV fc ) j 


t=i 


3 


EE 3 m - 1 (eW = 

k= 1 j=l i=l fc=l 


e'P'T = ^2^2 3 m Vd'T = 


^ 1 - pJ 

.7=1 p 


Note that moving the derivative inside of the summation is justified by Theorem 7.17 in 
[25] . at least for t in a closed interval [—5, 5] for small enough 5 > 0. Setting t — 0 gives us 

©. ’ □ 

Finding the moments of a random variable through its cumulants is a well-known tech¬ 
nique, but we show the details here for completeness. Note that Mx p {t ) = e 9p ^ is the 
exponential function composed with another function. Applying Faa di Bruno’s formula (J7J) 
to Mx p (t ) shows that 


k =1 

Setting t = 0, we arrive at the following result. 

Theorem 3.3. Using the notation from SectionUl and Lemma UTR for each p e (0,1) and 
for all me N, 

E(Xp = B m { K p, K ,f\--. ,4"*)). (15) 

The exact calculation of the moments in Theorem 13.31 can be quite unwieldy for large 
values of m. Thus, it’s often more insightful to observe simpler asymptotic formulas for the 
moments instead. In m the following asymptotics as p f 1 were proven for the input series 
of the Bell polynomials in 13.31 

p k 11 ln(l — p) 

1 — pr 1 — p 1 — p In p 

and, for j > 1, 

y- k j p k ~ jK(j + 1) _ jK(j + 1) 

1 - p k ~ (1 - p)i +1 ~ - ln J+1 p ’ 



(16) 

(17) 
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where Q(j ) = Yl'kLi 1 jW is the Riemann zeta function. 

As noted in the previous discussion on Bell polynomials, for each fixed m 6 N, the poly¬ 
nomial B m (x i,..., x m ) includes only one term of the form x™, and this term has coefficient 
1. Recalling the construction of the Bell polynomials (JH]) and the pigeonhole principle, it is 
clear all other terms of B m [x i,..., x m ) is a multiple of some variable other than X\. 

Note that from (TT6|) . ( ln (i^ p ) )" t (^p 1 ' > )'" —> 1 as p f 1. Also, it follows from that for 
all j e {2,3,4,... ,m}, 

ln p y K e) c? - i) ! cc?) 0 

ln(l-pp f -ln'(l-p) 

as p f 1. Since every term in (TT5|) is a multiple of some cumulant other than k p , except for 
the (ftp 1 '*)” 1 term, multiplying E(X™) by ( i n (mp) ) an d taking the limit as p f 1 eliminates 
all terms except for the afore mentioned term, where the limit converges to 1. Therefore, we 
can obtain the following asymptotic result: 

Corollary 3.4. For all m e N, E{X™) ~ (^^) m as p t 1. 

To further clarify the special role of the number Z p = 111 in our model consider 

Up := number of frogs who reached the site — Z p , 


where, for simplicity and clarity of the subsequent computation, we treat Z p as an integer 
number. Then, using (JHJ) and the Markov property, observe that 


E(U p ) = P(S n = — Z p for some n > 11 S 0 = x) = 


P 




x=0 


x=0 


p Zp 

i - p 


= l. 


This result can be heuristically interpreted as an illustration of the fact that — Z p serves as 
the most distant place, though barely, is still accessible to the frog population (in average 
only one frog can reach that far). 

The next corollary is merely a couple of special cases of Theorem 13.31 but will be useful 
in the discussion of the asymptotics of X p . 


Corollary 3.5. For any p e (0,1), the expected value of X p is given by: 


__ n x 

X=1 


Vy( 1) + ln(l - P) 

In p 


(18) 


where ip p is the q-digamma function defined by ([5]). Also, the variance of X p is: 


VarPQ 


V-'(I) 


X=1 


1 — p x In p 


(19) 


We now seek to determine the asymptotic of the distribution of X p when p f 1. Theo¬ 
rem 13.11 and the result in dT8|) both hint that X p grows at roughly the same rate as 
as p rises close to 1. The nature of this growth can be revealed by the asymptotic of the 
variance, which grows sufficiently slow to guarantee scaling limits for X p . 

















Theorem 3.6. Let Y p share the same distribution as -^~^X p for each p G (0,1). Then , as 
p t 1, Y p ► 1 in probability. That is, for all e > 0, 

lim P(|Yp — lj > e) = 0. 
p->i~ 

Proof. By (jTHD and (fTBlh E(Y p ) = hi | 1 1 1 ^^ E(X P ) —» 1. To show the convergence in probability, 
it is enough to prove that Var(Y^) — y 0 as p f 1 . By (1T9l) and (11711 . Var(X p ) ~ ~ YhEJy 

Thus, Var(V p ) = j^yVar (X p ) ~ 0. Since E(Y p ) -+ 1 and Var(V p ) -G 0, 

convergence in probability follows. □ 


Aside from the probabilistic implications of Theorem 13.61 we can also use the generating 
functions of Y p to construct some limit identities involving the g-Pochhammer symbol. Since 
Y p — y 1 in probability (and thus in distribution), the moment generating function of Y p con¬ 
verges pointwise to that of the degenerate variable at 1. Similarly, the characteristic function 
also converges pointwise to the characteristic function of the same degenerate variable. This 
observation leads to the following corollary: 


Corollary 3.7. For any z > 0, 

lim-- 

Ptl p- p) ^ 


Also, for any tel, 


lim ■ 

pti 


(p; p) c 


e^p;p) ( 


= e 


We wish to prove a stronger convergence of the sequence of random variables {Pp} pG ( 0 i). 
Theorem 13.61 implies that an appropriate discretization {V p71 } n eN of {Y p } pe ( 0il ) can be chosen 
to achieve the almost sure convergence. The following result identifies a class of sequences 
{PnjneN that ensures the almost sure convergence of the discrete sequence Y pn . 

Proposition 3.8. Let {p n }n &n Q (0,1) be a sequence such that p n f 1 and hi(1 ^ y G d 2 . Let 
{^pnlneN be a sequence of random variables in the probability space (fl,E,P) such that Y Pn 
has the same distribution as -^f^-^X Pn for each n. Then, as n —> oo, Y Pn —» 1 a. s. 

Proof. Let e > 0 be given. By the Borel-Cantelli Lemma, a sufficient condition for a. s. 
convergence is that P{\Y Pn — 1| > e) G LA. Choose p G (0,1) such that | E{Y p ) — 1| < |. 
Then, using Chebyshev’s Inequality, 


P(\Y p -l\>e)<P(\Y f - E(Y„)\ + | E(Y„) - 1| > e) 

< P (lU - E(Y„)\ > € - or | E(Y P ) - 1| > |) 

= p(|n-B(n)|>l)<ivar(y„). 

1 2 2 1 

As noted above, Var (Y p ) = Var( X p ) ~ ■ Replacing p with the terms of 

(PnjneN described in the theorem statement, we see that P(\Y pn — 1| > e) G i 1 , proving our 
result. □ 
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The class of sequences defined in the hypothesis of the proposition above includes those 
of the form p n = 1 — e ~ nC , where c > ^ is constant. For further research, we wish to broaden 
the class of sequences that lead to the a. s. convergence. For instance, instead of assuming 
that the models corresponding to different values of p are independent, one can consider a 
standard hierarchial coupling of the underlying random walks leading to the setting where 
P(X p i < X P2 ) = 1 if pi < p 2 . In that case one can consider for example p n = 1 — n~ c , c > |, 
and imitating Etemadi’s proof of the law of large numbers (see, for instance, Section 2.4 in 
m, namely first considering subsequences k(n) = [a n J with an arbitrary a > 1 and then 
using the fact that the ratio of lu | 1 ] 1 ^ ( ” ) ^ and ) converges to a c when n —> oo, prove 

the almost sure convergence of Y pn for the sequence p n = 1 — n~ c by finally taking a to 1. 


4 More general frog distributions 77 

Now let’s consider the frog model with drift in which its frog distribution r] = {r] x }^L 0 is a 
sequence of natural numbers, i.e., p x > 1 for all x > 0 and rj x — 0 elsewhere. Main results of 
this section are stated in Theorems 14.11 and 14.21 below. 

According to Theorem 2.1 in G3. in order to have transience in the frog model with 
drift p, and hence an almost surely finite minimum of its range, we must assume that 
Yl^=o r lxP x < oo. Since we will be dealing with a continuum of choices for p, it will be useful 
to designate all of the frog distributions that will guarantee transience in the frog model 
with drift. Hence, we define the following set of integer sequences: 

OO 

HI := 1 77 G pj Nu l°l : ^ r\ x p x < 00 for all p G (0,1) j. 

x=0 

It’s worth noting that HI does contain unbounded elements, such as 77 = {1, 2, 3, 4,... }. I 11 
fact, any integer sequence 77 such that p x = o(a x ) as x —> 00 for any a > 1 is in HI. 

Similarly to the single-frog case, we can construct a family of independent random vari¬ 
ables {X PtV ) pe f 0t i) tVCM on the probability space (H,^, P) such that for each p e (0,1) and 
rj G HI, X pv shares the same distribution as the negative of the minimum of the frog model 
with drift p and frog distribution 77 . 

By using similar ideas as in the single-frog case, we can find the distribution of X p ^\ 


p(x p>ri < x) = J](i - p x + k +y\ (20) 

k=0 


For simplicity, we will define the integer sequence {A fc }^T 0 by A 0 = po an d for all k > 1, 
Xk = Pk — Pk-i- The density of X p ri is then 


p(x m = X ) = n( i -p I+ ‘ +1 )”‘ 

k =0 


~[(i ~p x+k y k 

k =0 


00 00 

= U(i - p^+y* (1 - jj( 1 - P x+k ) A ^ . 

k =0 fc=0 


( 21 ) 
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One special case to consider is r] x = n G N for all x > 0. Then, (l20]i and (|2TT) become 


p(x„ < x ) = na - p x+t+, r = (/ +1 ;p)«, 

k =0 


P(X p „ = x) = (l-(l-pT)(P x+ 1 -,P)" x - 

With frog distributions 77 that differ from the single-frog case, we could assume that the 
moments of X Pt7J grow at different rates than hl | 1 ~ p ^ found in Corollary 13.41 However, by 
the theorem below, if 77 grows at a “slow enough” rate, the moments of X P)Tj will behave 
asymptotically similar to those of the single-frog case. 

Theorem 4.1. For each p G (0,1), let Z p = ln ^~ p ^ . Suppose {r]k}kLo £ H such that 
lim p | 1 (l — p) 1+s ET=oVkp k = 0 for all 5 > 0. Then, for all m, n G N, E(X™ V ) ~ Z™ as pf 1. 

The proof of the theorem is given below in this section, after a shirt discussion of the 
result. Note that according to the theorem, any even frog distribution rj = {n, n,.. . }, where 
n G N, will produce the same asymptotic rate for the moments. Not only that, but there 
exist unbounded choices for 77 that produce the same rate as well. One simple example of 
such a choice is 77 = {1, 2, 3,... }. 

A consequence of Theorem 14. II is an analogue to Theorem 13.61 in the single-frog case, which 
reveals that lll ^ i ~ p ' > is also an appropriate scaling limit for X p ^s convergence in probability 
as p 11. 


Theorem 4.2. Let 77 G HI satisfy the conditions in the hypotheses of Theorem \f.l 
share the same distribution as 
probability. 


ln P X 

ln(l-p)' V P,»7 


for each p G (0,1). Then, as p f 1, Y pr 


Let Y p p 
—> 1 in 


To clarify the intuition behind this result it is instructive to consider the case of an even 
frog configuration rj x = m for all x > 0, where m G N is a fixed integer, and observe that 
the corresponding model can be thought of as a composition of m independent models with 
rj x — 1. In this case, Theorem 14. 2 1 is a direct implication of the result in Theorem 13.61 following 
by a simple observation that since the random variable Y p is asymptotic to a constant, the 
same is true for its analogue Y ppi in Theorem 14.21 which is the minimum of m independent 
copies of Y p . From this perspective, Theorems 14.11 and Theorem 14.21 can be viewed as an 
indirect extension of this argument to sequences r] x growing sufficiently slowly, so that they 
can be well enough approximated by initial configurations with an even distribution of frogs 
(notice that the further is the initial placement of a frog from the origin the less relevant it 
is for the asymptotic of Y PtV ). 

Before we begin the proof of Theorem 14.11 we must first introduce a lemma. 

Lemma 4.3. For each p G (0,1), let Z p = Then, for all S > 0 and m gN, 


£(Z,(l+<S)+x)V~ 

X=1 


z'p 1 + S) 
1 - p 


(22) 


as p f 1. 


11 


















4-3. For the quotient of the left- and right-hand sides in (1221) we have 


Proof of Lemma 


d-p)E( 1 + 


X=1 


X 


oo m 


Z p ( 1 + 5) 


P x = 


a-dE^EI™) 

x=l j =0 7 


ad 


Z J p (l + 5)i 


(1 ^zia+syyj, 


m 


“ (1 p ^z}(i+sy\jJo-t>y +i 


E 

tc=l 

m \ j 


x°p x 


E 


m! 


^ (m-j)\ [Z p (l + 5)(1 -p)p" 


Now, as p f 1, —>■ —1 and ln(l — p) —* —oo. Thus, Z p (l — p) —>■ oo. Hence, for any 

j > 0 , the j th term of the above sum goes to 0 as p f 1. So when taking the limit, the only 
term in the sum that survives would be the 0 th , which is equal to 1 for all p. □ 

We now proceed with the proof of the above theorem. 

Proof of Theorem \4-l\ Since the frog model corresponding to X pv contains more frogs than 
the single-frog case, by Corollary 13.41 


lim inf 

pfi 


E(X™) 

Z™ 


> 1 . 


For the other inequality, note that for any p G HI, 


E(X™) < E(X%), 

where Ok — naa x{r)j : j = 1,2 ,,k}. Thus, we can assume WLOG that p is a nondecreasing 
sequence, and hence > 0 for all k. 

Choosing 5 > 0, we consider the sum for E(X™^) and split it at the point [Z p (l + 5)J. 
For the tail sum, we find that 

OO OO OO 

x m P(X p ,r, = x)< x m (l- JJ(1 - p x+k ) A 

x= L^p(lH-5)J +1 x= \_Zk =0 

oo oo 

~ E( z -(i +b m (i - IK 1 - /> z '< 1 + v+‘) a ‘). 

x=l k =0 

Since p Zp = (1 — p), we can see, after expanding the infinite product up to the first-order 
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terms, that 


x) <Y.( Z ,(1 + i) + x) m (l - n<l - (1 - rt 1+ V + ‘) A ‘) 

x=l k=0 

oo oo 

< E( z -(i+ j )+h m (E - p) i+ y +l ) 

X=1 k=0 

oo oo 

= (1 - p) 1+i (E VjEW 1 + i) +*)V 

/c=0 x=l 

oo 

~(i-rt 1+J (E’'‘^) z r( 1+i r- 

fc =0 

The equivalence result in last line comes from Lemma 14.31 combined with the fact that 

Eto^p k = vo + (i-p)EZiVkP k - 

We also derive the following upper bound for the finite sum: 

LZp (1+,5)J 

Y, xmp ( X w = x )< [Z p \ m {l + S) m P(X p , v < Z p (l + 5)) 

x=0 

< Z™(1 + 5) m . 

Thus, we can derive an upper bound for the following limit: 

jp( Vxri \ 00 

lim sup < lim(l + (1 + (1* p) 1+S ■ ^ 

Pti pfi V U ’ 

= (l + 5) m 

rpf v'm \ 

Since 5 > 0 is arbitrary, lim sup^ z £ v < 1, and this completes the proof. □ 

5 Initial configuration 77 supported on the whole Z 

In this section we will provide asymptotic bounds for the minimum of the frog model’s range 
under the assumption that 77 is supported on all of Z. The main result of this section is 
stated in Theorem 15.11 below. 

Up until now, we have assumed that there were no sleeping frogs on any of the negative 
sites. With this assumption, all of the frogs on Z + would eventually wake w.p.l, and we 
only needed to observe the collective minima of those frogs. However, when we consider 
the transient frog model with configuration r) supported on all of Z, we now have a random 
number of active frogs originating from the negative sites that have the potential of expanding 
the range. We begin to explore the moments of the minimum of this case by the groundwork 
laid in the previous sections for the 77 supported only on nonnegative sites. 

Fix any 71 e N. We will assume throughout this section that r/ x = n for any x G Z, that 
is exactly n frogs are initially placed at each site of Z. Our proofs in this section rely on 


U x™P(X„ = 

x= Lz p (i+5)j+i 
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the following description of the “avalanche structure” of the model. We refer to the frogs 
initially located in the nonnegative sites of Z as the “first wave”. If we just observe the 
nonnegative frogs, we can locate the left-most site visited by the frogs from the first wave. 
We consider the frogs on the negative sites down to the minimum ever visited by the first 
wave to be the “second wave” of frogs being activated. Tracking the left-most site visited 
by the frogs from the second wave, we designate a “third wave” of frogs activated between 
the subsequent minimums. We will continue to label these activated frogs in terms of waves. 
Since we assume a transient model, there will eventually be a final wave of frogs w.p.l that 
never venture any more to the left than their initial locations. In this section, let the negative 
of the left-most site visited by any of the active frogs be X p>n . 

The following theorem provides upper and lower bounds for the m th moment of X PtTl for 
any given m G N. While the bounds contain the familiar Z p = term from Sections [3] 

and SI they are not immediately obvious from the previous results. 

Theorem 5.1. The following holds for any m G N : 

(a) The m th moment of X p%n is bounded above by a function : (0,1) —> [0, oo) such that 

^ (23) 

(b) For any function 5 : (0,1) — > (0,1) such that lim^i S(p) = 0 and lim p ^i(l — p) 6 ^ = 0, 
the m th moment of X p ^ n is bounded below by a function ifs '■ ( 0 , 1 ) —> [ 0 , oo) such that 

Mp) ~ m! z 7 ex p{ ^ | as p tl. (24) 

Remark 5.2. An example of a function 6 that satisfies the conditions of Theorem I5.il (b) 
is 5(p) = {ln(jd—)} a for a constant a G (0,1) and p > 1 — e~ l (as long as we are 
only interested in the asymptotic as p f 1 , the values of 5(p) can be assigned arbitrarily 
for p < 1 — e~ x ). In this case, ifs(p) ~ m! Z™ exp{mnexp(| ln(l — p)| 1_Q )} as p f 1. 
For the sake of comparison with the upper bound, notice that the latter can be written as 

n 

4>{p) ~ Z™(-^) 2 exp{mnexp(| ln(l — p)|)} as pf 1 , and that the parameter a e ( 0 , 1 ) can 
be chosen arbitrarily close to zero. 

To motivate and clarify the intuition behind the coupling construction employed in the 
proof of Theorem 15.11 given below, we precede the proof by the following observation. 

Remark 5.3. Consider the model described in Section 0 namely rj x — 1 for x > 0 and 
r) x = 0 for x < 0. Fix any 5 > 0 and consider the probability P p j that no one of the frogs 
initially placed to the right of Z p ( 1 + 5) will ever reach zero. For simplicity and without loss 
of generality we will treat Z p { 1 + 5) as an integer. Then, since p z A l + & ) = (i _ p) 1+5 , 

OO OO 

= IE 1 - p Zp[l+5)+j ) = IE 1 - (! - P) 1+ V). 

3 =1 3 = 1 
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Since 1 — x > e 2x for all x > 0 small enough, we obtain that for p sufficiently close to 1, 


P p ,s >exp{-2(l-p) 1+5 ^^} = exp{-2p(l-p) 5 }, 

3 = 1 

and, hence 


lim P pS = 1. 

pti 


(25) 


Now, recall that the results of Section^ indicate tight concentration of the distribution of X p 
around Z p as pf 1. On the other hand, heuristically, (125]) indicates that for large values of 
p only the first Z p frogs are relevant to the dynamics of the model. To further support this 
claim, consider the probability Q P) & that no one of the Z p ( 1 — 5) frogs initially placed at the 
first Z p (l — 5) nonnegative integers will ever reach —Z P {1 — 5). Similarly as before, we treat 
Z p ( 1 — 5) as an integer. Then 


Q P ,s = n (l-p z ' <I - <w )<exp{-(l-p) 1 - ! Y, 

3=0 3=0 

= exp{—(1 - p)- s (l - = exp|—(1 - p)-‘( 1 - (1 - p) 1 - 5 )} 


< 


exp j —(1 - p)“ 5 (l - e - p ^- s) ) }. 


Thus 


lim Qps — 0. (26) 

pti 

Heuristically, (125]) along with (125]) tell us that the dynamics of the model considered in 
Section El is for large values of p similar to the dynamics of a modification where rj x equals 
1 only if 0 < x < Z p and is 0 otherwise. The proof of Theorem 15. II given below is using 
an interpretation of the model considered in this section as an “avalanche” of the models 
described in Section [7] and is using the above heuristic observation to produce upper and 
lower bounds of Theorem \5.1\ for the model range. 

We now turn to a formal proof of Theorem 15.11 

Proof of Theorem 15.11 

(a) We first provide an upper bound for the moments of X P)n by coupling the frog model 
with the following variant. First observe the minimum location reached by the frogs that are 
initially placed on the nonnegative sites and obtain a second wave of active frogs. Modify the 
original second wave in the following way. Put n more frogs on each site to the right of the 
minimum reached by the first wave and suppose that only the second wave can activate them. 
Note that without the consideration of activation times, this set of frogs is a translation of 
the configuration of nonnegative frogs considered in Sections [3] and El Find the minimum 
bound for this modified configuration and, for the resulting next wave, add frogs to all of the 
right-side sites in a similar fashion. Since there is a positive probability for the nonnegative 
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frogs never reach — 1 , and each wave of frogs is a translation of this case, the waves will 
terminate w.p.l. Let W p>n be the negative of the minimum bound for this model. 

A formal definition of the distribution of W pn can be given in the following manner. 
Define the sequence 77 = {rj x } xeZ G Z+ as follows: 

~ [ n if £ > 0 

y 0 otherwise. 

Such a configuration has been considered in Section [4] Let X PiU be the negative of the 

(k) 

range of the corresponding model and let k E Z, be independent copies of this random 
variable. Let 


T p , n = inf {k G N : Xj$ = 0} 


and 


W p ,n = 


V xw 

/ . p,ni 


k =1 


where, as usual, the empty sum (when T pn = 1) is interpreted as zero. 

To facilitate our computations we will actually use the following equivalent modification 
of this definition. Let 


£ P ,n = P{X P)n = 0) = ((p-p) oo)", 

where X p>n is the minimum of the range of the case with n frogs on each nonnegative site, 
introduced in Section [2 Thus e p = P (xjfil = 0) for any k e N. Let T p n be a geometric 
random variable with parameter e p . Namely, 

P(T Ptn = k) = (1 - £ p ) k £ p , k — 0,1,.... 


Notice that, limp_>.i£p in = 0, and hence T p>n converges to inhnity in probability as p t 1- 
Clearly, T p ^ n has the same distribution as T p ^ n — 1. Let y PtU = {Y p $} keZ be an i.i.d. sequence 
independent of T p ^ n and such that for any j G N, 


P(Y p , n - 3 )- P{X p , n - 3 \X p , n > 0) - 1 _ { p{ P x pn = 0 ) 

Finally, let 

T 

± Pi n 

vv P,n 1 p,n ' 

k= 1 

Clearly, W p . n has the same distribution as W p - n . 


16 



For m E N, we look at the m th moment of W pn conditioned on T pn : 


E(W™\T p , n ) = E[[Y,Z Y iS 


T, 


p,n 


< 


eZ' e ko”] 


mi 1/m 


where we use Minkowski inequality and the fact that Y pp } are independent of T p n . From this 
conditioned expectation, we approximate the m th moments of W PtU as p 1 1: 


E(W^ n )=E{E(W^\T nn ))<E[{T p , r ,) m }-E[(Ym) m ] 


(27) 


where Z p — as in Section QI To evaluate the moments of T pn we used the following 

known result whose short proof is supplied for reader’s convenience. 

Lemma 5.4. For e e (0,1), let T e be a geometric random variable with parameter e p . 
Namely, 

P(T e = k) = (l-£) k £ , A: = 0, 1, .... 

Then, for any m e N, E(T™) ~ £~ m as £ f 1. 

Proof of Lemma fFf \ We have as e f 1, 

°° /»oo pc 

E(T T) = ^P(Tf >*0= / P(Tr>x)dx= / 

b .—n 0 J 0 


P(T e >y)-my m - L dy 


k =0 
eoo 


= m( 1 — e) 


(l - £) y+l ■ my m ~ l dy = m(l - £) / dy 

Jo 

r (m) 


m! 


ln(l-£)| m 
The proof of the lemma is complete. 


□ 


Using (127)) . we finally arrive at the asymptotic bound in (a) through the asymptotic of 
the g-Pochhammer symbol derived from [[14]: for p = e _t , as t 0, 


(p; P)c 


2vr 



exp 


- 7 T 


6 1 


Note that t = — lnp ~ 1 — p, and we have (a). 


□ 


(b) Now, consider another variant of the frog model. Define a function 8 (p) : (0,1) —> (0,1). 
To simplify notation, we will occasionally use 8 p for 8 (p). For a given p, consider the config¬ 
uration fj with rjk = n if k £ {0,1, 2,..., Z p ( 1 — 8 P )} and rfr. = 0 elsewhere. For simplicity 
and without loss of generality, we will assume that Z p (l — 5 P ) is integer-valued. Start the 
model, and see if the frogs eventually reach the site —Z p (l — 5 P ). If they do, activate n frogs 
on each of the Z p (l — 5 P ) sites to the left of the origin. Now observe if the newly activated 
frogs reach the site —2Z P (1 — S p ). If they do, activate n frogs on all the Z p (l — S p ) sites to 
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the left of the sites previously activated. Continue this procedure indefinitely. Note that at 
any observed time step, this model will always have fewer active frogs than the frog model 
with initial configuration 77 with 77 *, = n for all k 6 Z. By Theorem 2.1 in j!2j . the frog model 
with this configuration is transient. Hence, the variant model is transient, and the process 
will eventually stop producing new sets of frogs from the left. Let V Pj s be the negative of the 
minimum of the range of the variant, and let r p _s be the number of activated blocks of the 
length Z p (l — 8 P ), not including the initial one located within Z + . 

Let 

=''“n ‘(i - pk (1 

7=0 


By this definition, 6 P) s is the probability that none of the frogs located at the first Z p ( 1 — 5 P ) 
nonnegative sites will ever reach the half-line on the left of Z p (l — Sp). The event that a given 
block of frogs nZ p (l — d p ) active frogs reachs its goal and produces another active frogs block 
of length Z p ( 1 — 5 P ) is independent of that of any other set. Hence, viewing the event of a set 
reaching its destination Z p (l — S p ) to the left as a “failure”, we have that t P: s is geometrically 
distributed with parameter 9 p j. Namely, 

P(t p ,s = k) = (1- 9 P! s)%,s, k = 0 , 1 ,.... 

Notice that since p Zp = 1 — p, 


Z p (l—5p)—l 


6 p,5 = 


Zp(l—5 p ) — l 


= exp 1 —n 


n (1 - (1 - < exp(-r»(l - p) 1 - 5 ' Y, y 

7=0 7=0 

l — P Z p(. 1 ~ 6 p)\ / 1 — (1 — p) 1-<5 ' 


P 

(i - P y 


= expl — n 


(1 -pY 


With the 6 chosen with the constraints specified in the theorem’s hypotheses, we see that 
9 p $ — y 0 as p t 1. Hence, by Lemma 15.41 E(t™ s ) ~ m\9 as p t 1- Clearly, X pn is 
stochastically dominated from below by V Pt ft = T Pt s ■ Z p (l — S p ). The lower bound for the 


moments of X PtU is therefore 


yyyy(l-t)" ~ r ^-Z™~m\Z™ ew {mn l - 

~ ml Z™ exp|- 

as p t 1. The proof of the theorem is complete. 


1 — s 0 


mn 


(i - py 


□ 
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